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Non-resonant diagrams in radiative four-fermion processes 

J. Fujimoto, T. Ishikawa, S. Kawabata, Y. Kurihara, Y. Shimizu '^ and D. Perret-Gallix^ 
'^ Minami-Tateya Collaboration, KEK, Japan 
'=LAPP-IN2P3/CNRS, France 

The complete tree level cross section for e^e~ -^ e^UeUdj is computed and discussed in comparison with 
the cross sections for e^e~ — > e^ Ugud and e^e~ — > udud. Event generators based on the GRACE package for 
the non-radiative and radiative case are presented. Special interest is brought to the effect of the non-resonant 
diagrams overlooked so far in other studies. Their contribution to the total cross section is presented for the 
LEP II energy range and for future linear colliders {^/s =500 GeV). Effects, at the W pair threshold, of order 3% 
[e'VeUd) and 27% (udud) are reported. Similar behaviour for the radiative case is shown. At ^/s — 500 GeV, the 
relative contribution of the non-resonant diagrams for the radiative channel reaches 42.5%. 



1. Introduction 

Two years ago, a complete calculation of 
two four-fermion final states in e+e" collisions, 
e+e" -^ e~Ueud and e+e" -^ udud, was re- 
ported at the Sotchi meeting ||l[. It was shown 
that the non-resonant diagrams play an impor- 
tant role below the W pair threshold (up to 27% 
at ^/s = 150 GeV) and exhibit a non-negligible 
effect at higher energy (7% at i/s = 190 GeV). 
In these results, in addition to a set of exper- 
imental loose constraints, a cut on the 9 angle 
of the outgoing electron was applied to cope with 
gauge violations appearing in the subset of 7 — 14^ 
diagrams. Below threshold, the non-resonant di- 
agrams competing only with the off-shellness of 
the resonant ones, their relative contribution be- 
comes quite large. The non-resonant diagrams 
are dominated by t-channel j — W graphs in sim- 
ilarity with the well-known two-photon (7 — 7) 
processes. Hence the cross section increases with 
the energy and the contribution of non-resonant 
processes keeps growing in contrast to that of res- 
onant diagrams which falls with the energy. This 
effect was confirmed at this workshop by the anal- 
ysis of selected final states based on a dedicated 
four-fermion generator g . The non- resonant pro- 
cesses are essentially single W or no W produc- 
tion. Some examples of these type of diagrams 
for the radiative case can be seen on Fig. lb, 2. As 



far as the W mass reconstruction is concerned, 
these diagrams should be counted as background 
processes, and they will contribute to the overall 
mass and width measurement uncertainties. 

It has been recognized, since some time [pj, 
that the radiative corrections in W pair produc- 
tion have to be included for a precise comparison 
with measurements at the e'^e~ colliders. Initial 
state radiations (ISR) reduce the available center 
of mass energy. Close to the W pair threshold, 
ISR will, therefore, enhance the relative contri- 
bution of non-resonant diagrams. Final state ra- 
diative corrections will contribute to additional 
uncertainties on the mass reconstruction, unless 
the photon is detected and properly accounted 
for in the mass reconstruction algorithm. How- 
ever, photons coming from the intermediate W, 
which can be mistakingly recognized as final state 
bremsstrahlung, will introduce additional uncer- 
tainties. 

Although a number of studies has been car- 
ried out in these last years, including 0[g, we 
present, for the first time, the complete tree level 
computation of one of the radiative four-fermion 
processes. This work can be used as a benchmark 
for testing the validity of simplified and therefore 
faster generators. 
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Figure 1. Some examples of diagrams from the 
radiative process e+e^ — > e^Deud'j: a) Resonant 
diagrams, b)Non-resonant j ~ W diagrams. 

2. The GRACE system 

The GRACE system Q has been used to 
perform the very lengthy computations involved 
in the study of most radiative processes. The 
GRACE package is a complete set of tools for 
computing tree level processes. All the phases 
involved in a given computation are covered: 
from the process definition (specification of the 
initial and final particles) to the event genera- 
tor. It is composed essentially of three compo- 
nents: the diagram generator, the construction of 
the matrix element based on helicity amplitude 
function from the CHANEL library and the 
multi-dimensional phase space integration pack- 
age BASES Pi associated with the event genera- 
tor SPRING g. 

Fermion masses are properly introduced in the 
helicity amplitudes. The particle width is intro- 
duced into the gauge boson propagators when 





Figure 2. (Fig.l cont.) Some examples of 
diagrams from the radiative process e+e^ —f 
e~'Deud'y: Other non-resonant diagrams. 

the denominator can vanishes for positive squared 
momentum transfer. A gauge invariance checking 
program is automatically built by the system. 

Only the kinematics for complex processes has 
to be written by the user. The issue, here, is to 
figure out what is the best transform to apply to 
the basic integration variables to regularize the 
divergencies appearing in the integration of the 
matrix element. The best kinematics will give 
the best accuracy for a given computation time. 

3. Four-fermion production 

The results presented hereafter have been ob- 
tained using the following set of parameters: 



Mz -~ 


= 91.1 GeV 


Tz -~ 
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= rud = 0.1 GeV. 



The width of the W is taken from the Particle 
Data Table; Tw = 2.25 GeV. The gauge boson 
(W,Z) widths are assumed to be constant in the 
calculation. Furthermore some realistic experi- 
mental cuts have been introduced: 
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where 9e,0-y and 9q^q are angles measured from 
the incident e~ beam and E^^q^q are the energies 
of final e^,(7 and q, respectively. The electron 



polar angle cut is introduced to avoid an undesir- 
able gauge violation for extremely forward direc- 
tions y] caused inevitably by the finite width of 
W boson when the subset oi j — W diagrams is 
considered. 

The gauge invariance of the amplitude is 
checked numerically for a random selection of the 
boson gauge parameters at several points in the 
phase space. The errors are within the preci- 
sion of numerical calculation (typically less than 
0(10^^^) in double precision) when the W width 
is turned off. Further gauge cancellations are 
checked to confirm that the obtained cross sec- 
tions are stable against finite width. The pro- 
cedure is as follows: first extract the product 
of gauge boson propagators without width as an 
overall factor of the whole amplitude. Then re- 
place them by those with finite width. As a result, 
one is left with rigorously gauge invariant ampli- 
tude[pOt. For the resonant diagrams, this method 
gives the same results as ordinally method. While 
for the ^ — W diagrams, this gives 10% smaller 
results. 

3.1. The non-radiative case 

The contribution of the non-resonant diagrams 
has been considered for two final states, e~v^ud 
and udud. Fig. 3 shows the cross sections ver- 
sus the center of mass energy around the W 
pair threshold. The solid line represents the case 
where all diagrams are included, the dashed line 
the contribution of the resonant diagrams only. 
In the udud channel, one can observe the onset 
of the Z resonant diagrams as a small shoulder 
starting around 2Mz. The behaviour below the 
W pair threshold can be better seen in fig. 4 where 
the ratio of the three main resonant diagrams over 
the complete calculation is plotted for the same 
range of center of mass energy. In the e~Deud 
channel, the effect is of the order of ±3% around 
\/s = 160 GeV, while in udud it reaches —28%. 
At the energy of future LC, 500 GeV, the contri- 
bution of non-resonant diagrams is quite signifi- 
cant in e^Veud but is negligible in udud. 

Similarly to the 7 — 7 events, the non-resonant 
processes are predominantly peaked in the for- 
ward/backward direction. Fig. 5 shows the be- 
haviour of the cross section versus the electron 



0. 




180 200 
VS (GeV) 



Figure 3. Cross sections for a) e+e~ — » e~Ueud 
and b) e^e" -^ udud versus CM. energy with 
9e > 8°. The solid curves are for all diagrams 
and dashed ones for three resonant graphs. The 
dotted line in a) is the contribution from the 7 — 
W and Z — W processes. The estimates at 500 
GeV are indicated by arrows. 



polar angle cut. The non-resonant contribution 
vanishes as the cut is raised. 

3.2. The radiative case 

In the covariant gauge, the GRACE package 
generates 752 diagrams when all the interactions 
are included, but actually 142 diagrams are used 
in the unitary gauge with Higgs-fermion coupling 
being set to 0. The full tree level computation 
requires, for a decent precision, some 100 hours 
on a HP 735 computer. Thanks to the vectoriza- 
tion of the integration package BASES, the most 
extensive computation could be performed on the 
KEK vector processor within 1 hour. Fig. 6 shows 
the contribution of the non-resonant terms in the 
total cross section versus the CM. energy. The 
behaviour is similar to the non-radiative case. In 
fig. 7, one can see more clearly the contribution 
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Figure 4. Energy dependence of the relative con- 
tribution of resonant diagrams in a) e+e^ -^ 
e~Deud and b) e+e^ -^ udud with O^ > 8° to 
the fuU computation. 

of non-resonant diagrams similar to fig. 4 for non- 
radiative case. 

The dependence of the cross section on the elec- 
tron polar angle cut 9e is presented in fig. 8. Fig. 9 
shows the distribution of the invariant mass of the 
pairs of final state fermions. The general shape is, 
as expected, a Breit-Wigner distribution. How- 
ever some background coming from the onset of 
non-resonant processes can be seen. In addition, 
for the radiative process (fig. 10), a mass shift of 
the order of 1 GeV is observed on the e~De in- 
variant mass due to final state radiation. 

4. Conclusions 

A complete calculation of the radiative four- 
fermion final state e"'"e~ -^ e^Dgud'y has been 
performed, for the first time, using automatic 
Feynman diagram calculation techniques imple- 
mented in the GRACE system. Comparison with 
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Figure 5. Electron angular cut dependence for 
e+e- -^ e-Deud (CM. energy = 180 GeV). Solid 
line: all diagrams, dashed line: three resonant 
diagrams and dotted line: j — W and Z — W 
processes. 

two non-radiative computations e~i'eud and udud 
has been carried out. The contribution of the 
non-resonant processes is large below threshold 
(at least for udud channel) which makes the W 
mass determination from the excitation curve 
more difficult. At high energy, the contribu- 
tion increases, reaching 46% {e~Peud) and 42.5% 
(e^Deud'y) at y/s = 500 GeV. Determination of 
the W mass from direct reconstruction will have 
to deal with the intrinsic background from the 
non-resonant processes prominently in the for- 
ward direction. Furthermore the final state radia- 
tive correction will introduce a shift in the mass 
of the W unless the photon is properly taken into 
account in the algorithm of mass reconstruction. 
However, this last point may prove to be diffi- 
cult as the photon may be produced either from 
the initial state, the intermediate charged boson 
or the final state fermions. Only this last class 
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Figure 6. Total cross section for e+e^ -^ e^Peudj 
versus the CM. energy with 9e > 8°. The sohd 
hne represents the complete cross section while 
the dotted line the resonant term only. 
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Figure 7. Energy dependence of the contribution 



of non-resonant diagrams in 



^I'e'U-d'J 



with 6e > 8°. The ratio of contribution of reso- 
nant diagrams versus all diagrams is shown. 



plays a role in the W mass determination. More 
detailed studies are needed to check if additional 
kinematical cuts can improve the mass determi- 
nation. 

However to quantify the actual relative contri- 
bution of radiative events in the total cross sec- 
tion, soft photon and one-loop contributions must 
be taken into account. Large negative contribu- 
tions are expected, reducing, consequently, the 
number of radiative events and softening the mass 
reconstruction problem. 
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Figure 9. Invariant mass for non-radiative process 
in solid line for the ud system and in dotted line 
for the e^De system with 0e > 8°. 
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Figure 10. Invariant mass for radiative process in 
solid line for the ud system and in dotted line for 
the e^De system with 9e > S°. 



